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Abstract 

Analytic solutions of the Teukolsky equation in Kerr geometries are presented in the form 
of series of hypergeometric functions and Coulomb wave functions. Relations between 
these solutions are established. The solutions provide a very powerful method not only 
for examining the general properties of solutions and physical quantities when they are 
applied to, but also for numerical computations. The solutions are given in the expansion 
of a small parameter e = 2Mu, M being the mass of black hole, which corresponds to 
Post-Minkowski expansion by G and to post-Newtonian expansion when they are applied 
to the gravitational radiation from a particle in circular orbit around a black hole. It is 
expected that these solutions will become a powerful weapon to construct the theoretical 
template towards LIGO and VIRGO projects. 
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1 Introduction 



There are growing interests in analytic solutions of the Teukolsky equation Jl|] in the 
Schwarzshild and the Kerr geometries in the connection of the gravitational wave cos- 
mology. Since Teukolsky proposed the master equation for massless fields in the Kerr 
metrics, many efforts have been made to obtain the analytic solutions. The analytic ex- 
pressions valid for low frequencies were found by Page0], Starobinsky and Churilovj3| by 
matching the approximate solutions valid near horizon and far from it. Leaver |3J made 
the systematic study to obtain the analytic solutions of Teukolsky equation in the form 
of series of various functions. He found the solution in the form of series of Coulomb wave 
functions which is valid in the region far from the horizon and established the relation 
between that solution and the one in the form of the Jaffe type series which is valid near 
the horizon. 

Recently, Tagoshi and NakamuraQ determined numerically the coefficients of the 
post-Newtonian expansion of the gravitational radiation by a particle traveling a circular 
orbit around a Schwarzshild black hole. Sasaki proposed a method of post-Newtonian 
expansion to solve the homogeneous Regge- Wheeler equation by using Bessel functions. 
Subsequently, the extensive study on this line was made by Tagoshi and Sasaki and 
the result was compared with the one by Tagoshi and Nakamura. The application of this 
method to the Kerr geometries was made by Shibata, Sasaki, Tagoshi and Tanaka|J. 
Various other applications was discussed by Poisson and Sasaki |J. Now the problem to 
obtain the analytic solutions and the examination of their behaviors in low frequencies 
became the important and the urgent topic. 

In this paper, we report that we obtained the analytic solutions of Teukolsky equa- 
tion in Kerr geometries in the form of series of hypergeometric functions and Coulomb 
wave functions. The series solution of hypergeometric type is shown to be convergent 
in the region except infinity, while the one of Coulomb type is convergent in the region 
| x |> 1, where x = (r+ — r)/2Myl — (a/M) 2 with r+, M and a being the position 
of the outer horizon, the mass and the angular momentum of Kerr black hole. We es- 
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tablished the relation between two solutions with different regions of convergences. The 
solutions are interesting not only for the investigation of general properties of solutions 
as mathematical physics, but also are for various applications to the gravitational wave 
cosmology The solutions are essentially given in the e expansion which corresponds to 
the Post-Minkowskian G expansion and also corresponds to the post-Newtonian expan- 
sion when they are applied to the problem of the gravitational radiation from a particle 
in cicular orbit around a black hole so that our solutions are quite powerful to examine 
the e behaviors of various physical quantities. Our solutions are expected to become a 
powerful machine for numerical computation also because the convergences of series are 
well known. Thus the solutions will become a powerfull weapons for the construction of 
the theoretical template towards the gravitatinal wave observation by LIGO and VIRGO 
projects. 

Our work was motivated by Sasaki's work. We tried to improve his method for the 
solution of Regge- Wheeler equation because his method has several disadvantages: (1) it 
is difficult to obtain the higher order terms of e = 2Mu, where M and u being the Black 
hole mass and the angular frequency, (2) the expansion is not really the Bessel expansion 
because coefficients are also variable dependent and (3) the convergence of the series was 
unknown. In order to improve these difficulties, we considered the solution in the form 
of series of hypergeometric functions for the solutions of Regge- Wheeler equation and 
also for the Teukolsky equation in Schwarzshild metrics and showed that the coefficients 
of series can be determined systematically in the expansion of e due to the recurrence 



relations among hypergeometric functions which we found 1 10]. This solution is valid near 
the horizon and not at infinity so that away from the horizon we have to consider the 
solution in the form of series of Coulomb wave function which was found by Leaver Q. By 
matching these two solutions in the intermediate region, we obtained a good solution in 
the entire regions. After finishing our work, we happened to see the paper by OtchikfllT 



who discussed the analytic solutions of the Teukoksky equation in the form of series of 
hypergeometric functions and Coulomb wave functions. We found that our method is 
essentially identical to Otchik's method, but our solutions disagreed with his ones. We 
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compared our solutions and his ones and found that although various formulas which he 
presented were incorrect, his story itself turns out to be true. Since our results are all 
different from those by OtchikJTl[] and the results themselves are quite important for the 
application, we present all results in this paper. 

We start from the Teukolsky equation which is separated by writing 

V = e-**j m +S?(e)R uhn {r)- (1-1) 

The equation for R is 



AR" + 2(r — M)(s + 1)R' + 



K 2 — 2is(r — M)K 

h Aisuor - A 



R=Q, 



;i.2) 



where M is the mass of the black hole, aM its angular momentum, A = r 2 — 2Mr + a 2 = 
(r — r + )(r — r_) with r± = M ± V M 2 — a 2 where r + and r_ are positions of outer and 
inner horizons, respectively, K = (r 2 + a 2 )c<j — am, A = — s(s + 1) — Imauj + a 2 a; 2 . 
The function is the spin weighted spheroidal harmonics which will be discussed in 
Appendix A with the eigenvalue E ^2\ 

E = l(l + l)-^^ + [H(l + l)-H(l)-l\e + 0(e), (1.3) 

^ (/) - (2/-l)/3( 2 / + l) ' (L4) 
where £ = acu and I is the angular momentum which takes an integer or half-integer 
number which satisfies I > max(| m |, | s |). 

In Sec. 2, we give the discussion about how we arrive the analytic solutions in terms 
of hypergeometric functions and discuss their properties. In Sec. 3, the analytic solutions 
in terms of Coulomb wave functions are given following the work of Leaver. The relation 
between two solutions in two different convergence regions is established in Sec. 4. The 
low frequency expansion of these solutions is discussed in Sec. 5. In Sec. 6, the summary 
and remarks are given. 
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2 Analytic solution in the form of series of hyperge- 
ometric functions 



The radial Teukolsky equation has two regular singularities at r = r± and an irregular 
singularity at r = oo. In order to obtain the solution in the form of series of hypergeo met- 
ric functions, we have to deal with these regular singularities. Following the discussion 
in Appendix B, we take the form of R which satisfies the incoming boundary condition 
on the outer horizon. In particular, we choose the form given by (a_, (3+) in the notation 
in Appendix B with the variable x = u(r + — r)/eK as 

Rl n = e^i-xy^il-xY^pUx), (2.1) 

where e = 2Mtu, q = f ! ,K = ^l^q 1 and r = 
Then, the radial Teukolsky equation becomes 

x(l - x)pl" + [l-s-ie-ir-(2- 2it)x}pI' + ( v + lT ){v + I - ir)p v m 
= 2ien[-x(l - x)p u in ' + (l-s + ie- ir)xp v in ] 

+ [-A - s(s + 1) + v{v + 1) + e 2 - ie/c(l - 2s)]p u iw (2.2) 

Here we introduced the parameter v as the renormalized angular momentum which sat- 
isfies v — I + 0(e). Then, the right-hand side of Eq.(2.2) is of order e so that this form 
of equation is suitable to obtain the solution in the expansion of e. The zeroth order 
solution of Eq.(2.2) is the hypergeometric function. 

From the structure of the above equation, the solution may be written in the form of 
series of hypergeometric functions as 

oo 

tin(*) = E <Pr+v(x), (2-3) 
n=— oo 

where 

Pn+v(x) = F(n + v + 1 — it, — n — v — it; 1 — s — it — ir; x), (2.4) 

with the use of the renormalized angular momentum v rather than I. We expect that 
the series will coincide to the e expansion. In order for the coeficients of series (2.3) 
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to be a solved, it is essential that the coefficients a^'s satisfy the three term recurrence 
relation. For this, terms such as x(l - x)p' n+v and xp n+u must be expressed as linear 
combinations of p n+v+ i, Pn+u and p n+u -i. Amazingly enough, we found the following 
recurrence relations, 

(n + v + 1 — s — ie) (n + v + 1 — ir) 



•EPn+v 



+ - 



2(n + z/+l)(2n + 2z/ + l) 
ir(s + ie) 



-Pn+u+l 



(n + v){n + v + 1) 
(n + v + s + ie) (n + v + ir) 
2(n + v)(2n + 2is+ 1) 



(2.5) 



- X )Pn+u 



(n + z/ + ir) {n + z/ + 1 — ir) {n + u + 1 — s — ie) 
2(// • // • 1)(2// • 2// • 1) 
ir(l — ir) 



Pn+v+l 



+ 2^ S + k 



1 + 



Pn+u 



(n + v)(n + v + 1) 
(n + z/ + 1 — ir) (n + z/ + ir) {n + z/ + s + ie) 



Pn+v-li 



(2.6) 



2(n + z/)(2n + 2z/+ 1) 

which enables us to obtain the three term recurrence relation among a v n . The above 
recurrence relations among hypergoemetric functions can be proved by using the power 
series expansions. By substituting the form in Eq.(2.4) into the radial Teukolsky equation 
(2.2), we find that p v in becomes a solution if the following recurrence relation is satisfied: 



a 



n a n+l ' Pn a n ' 7n°n-l — U ) 



(2.7) 



where 



ien(n + v + 1 + s + ie) (n + v + 1 + s — ie) {n + u + 1 + ir) 
(n + z/+ l)(2n + 2z/ + 3) ' 



(2.8) 



Pn = -\-s(s + l) + (n + v)(n + is+l) + e 2 + e(e-mq) 
e(e — mq)(s 2 + e 2 ) 



+ 



(n + v){n + v + 1) 



(2.9) 



ien{n + v — s + ie)(n + v — s — ie)(n + v — ir) 
{n + z/)(2n + 2v - 1) ' 



(2.10) 
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By introducing the continued fractions 



Rniy) = -5 s ", L n {y) = (2.11) 



we find 



RM = — r^, LJv) = ^ — . (2.12) 

From these equations, we can evaluate the coefficients by taking the initial condition 
Oq = 1. The renormalized angular momentum v is determined by requiring that the 
coefficients obtained by using R n {y) agree with those by using L n (u), that is, by solving 
the transcnedental equation for v 

R n {y)L n ^{y) = 1. (2.13) 

If Eq.(2.13) is satisfied, we find 

i imri ^L = _ i im n <^ = ^_ / 214) 

n ^°° <-i n ^-°° <+i 2 
From the large n behavior of hypergeometric functions, we find by using the recurrence 
formula of hypergeometric functions (2.5) as [|l3f 

Jim Ws) = ]im Pn+u(x) = ! _ 2x + ((i _ 2x ) 2 - I) 1 ' 2 . (2.15) 

(x) p„ +i ,+i (x) 

From Eqs.(2.14) and (2.15), we find 

lim n < Pn +^ X ) = _ Jim KX + ,0r) _ 2x + ((1 _ 2x) 2 _ ^1/21 (2 lg) 

"- > - °°<+iPn+ v +i(a;) 2 

Thus the series converges in all the complex plane of x except for x = 00. 

As for the recurrence relation (2.8), we find aC^ _1 = 7^ and 7l^~ 1 = a v n so that 

aZ 1 ^ 1 satisfies the same recursion relation as that which a v n does. Thus if we choose 

a o = a o u l = 1) we have 

< = aX 1 - (2-17) 

Also, we find 

R n {-v - l)L n „ x {-p - 1) = R- n+l {p)L_ n {v) = 1, (2.18) 



which means that if v is the solution of Eq.(2.13), then — v — 1 is also the solution. 

It is easily seen that the solution R v in is symmetric under the exchange of v with —v — \ 
as follows. By using the formula 

„ (r) - r(l- 3 -»e-tr)r(2n + 2i/+l) +v+lT 
Pn+u{ } " T(n + u+l-tr)T(n + u + l-s-te) { ' 

1. 

x F(—n — v — z't, — n — v + s + ze; — 2n — 2z/; — ) 

x 

+ r(l - g - ze - zr)r(-2n -2u-l) x) - n -u + ir 
r(— n — z/ — zr)r(— n — v — s — it) 

x F(n + z/+ 1 -z'r,n + z/+ 1 + s + z'e; 2n + 2z/ + 2; -), (2.19) 



we can show 



iC = i?o + ^o , (2-20) 



where 



K = e-(-x)— £ fl , r(l- 3 -» e -ir)r(2n + 2 y + l) 

^ n r(n + z/+l-zr)r(n + z/+l-s-ze) 



n=— oo 



x (-x) ri F(-n-z/-zr, -n-z/ + s + ze;-2n-2z/;-). (2.21) 

x 

The behavior of i?^ n on the outer horizon (x = 0) is 

oo 

R- n - (-*r^ (e+T) E <C (2-22) 



n=— oo 



which gives the normalization of our solution. 

We can also show that Rq and Rq^ -1 are solutions which are independent each other. 
To see this explicitly, we consider the solution which satisfies the outgoing boundary 
condition on the outer horizon. The outgoing solution which corresponding to (a+,/3-) 
in the notation defined in Appendix B can be written as 

Km = e^(-x)^\l - xr s -^V uf (2-23) 
Now we expand p v out as 

oo 

P U out{*)= E <Pn+u(x), (2.24) 



where 

p n+v (x) = F{n + v + 1 + ir, —n — v + it; 1 + s + ie + it; x). (2.25) 

Similarly to the solution satisfying the incoming boundary condition, we find that the 
above series becomes a solution if the following recurrence relation is satisfied; 

«+i + P u n< + 7X-i = 0, (2.26) 

~u (n + v + l-s- ie)(n + v + 1 - if) . , 

d„ = -, rz rOt„, [1.2,1) 

(n + v + l + s + ie)(n + v + l + iT) n 

~v (n + v + s + ie)( y n + v + iT) 

In = 7 — : -r, - — ; ^7„, (2.28) 

(n + v — s — ie) [n + v — it) 

where a v n , (3^ and 7^ are defined in Eqs.(2.8)-(2.10). By inspection, we see that this recur- 
rence relation is reduced to the one in Eq.(2.7) by redefing systematically the coeficients 

as 

~ v _ Y{v + 1 - s - ie)T{v + 1 - ir)Y{n + v + 1 + s + ie)T(n + v + 1 + It) v 
fln ~ T(v+l + s + it)T(v+l + iT)T(n + u+l-s-it)T(n + u + l-iT) aw { ' } 

where we chose = 1. Now we take a zV ~ l = a$ = 1, then after some computation we 

find 

Rout — A u Rq + A^ u ^ 1 R I/ ~' 1 , (2.30) 

where 

_ r(l + s + ie + ir)T(u + l-s- it)T(u + 1 - it) 
v T(l-s-ie-iT)T(u+l + s + ie)T(u + l + tT)' [ ' ' 

This relation explicitly shows that Rq and R$ u ~ l are independent solutions of Eq.(2.1). 

3 Analytic solutions in the form of series of Coulomb 
wave functions 

Analytic solution in the form of series of Coulomb wave functions are given by 
Leaver 0. Here, we follow the discussion in Appendix B and start the parameteriza- 
tion to remove the singularity at r = r_. By using a variable z = u(r — r + ) = —ekx, we 
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Rc = z- 1 - s (l + -Y e ~ T Uz). (3.1) 



take the following form 

m. = r- 1 - 3 ^ 1 -u 

z J 

Then, we find 

z 2 fl + [z 2 + 2(e + is)z-u(u + l)]f v 

ffii,t\, ft , , • • \ fi eK(l + s + ie)(l-ir) 
= -eKz(f v +f v ) + eK(l + s + ie-iT)f,, /„ 

z 

+ [\ + s(s + I) - v(v + I) - 2e 2 + emq - en(e + is)]/,,. (3.2) 

If we consider v to be v = I + 0(e), the right-hand side of Eq.(3.2) is the quantity of 
order e so that this equation is a suitable one to obtain the solution in the expansion of 
e. 

Here we aim to obtain the exact solution by expanding f u (z) in terms of Coulomb 
functions with the renormalized angular momentum u, 

oo 

/„= £ b»F n+u (z), (3.3) 

n=— oo 

where F n+U is the unnormalized Coulomb wave function, 

F n+U = e-"{2z) n ^z T{n ^ n ± l - ~ + * + ie) <D(n + v + 1 - s + ie, 2n + 2z/ + 2; 2i*). (3.4) 

It is essential for the solution of Coulomb wave functions to be related to the one of 
hypergeometric functions, the renormalized angular momentum v takes the same value 
for both cases. 

By substituting Eq.(3.3) into Eq.(3.2) and using the recurrence relations satisfied by 
the Coulomb wave functions, 



1 (n + v+l + s — ie) is + e 

~z n+v ~ (n + z/ + l)(2n + 2z/+l) n+v+1 + (n + v)(n + u + 1) 



(n + v-s + ie) 
(n + v)(2n + 2u + l) 



{n + v){n + v + 1 + s-ie) is + e 

n+u (n + v + l)(2n + 2i/ + 1) n+!/+1 + (n + i/)(n + v + 1) n+ " 



(w + t/ + l)(n + z/-s + ie) 
(n + v)(2n + 2v + 1) 



we obtain the three term recursion relation of coefficients b^, 

<K + i + « + yXi = o, (3.7) 



, v . (n + v + 1 - s + ie) (n + v + 1 - s — ie) v 



"(n + v + l + s + it)(n + v+l + s-ie) an 

, v _ {v + n + s + ie){n + v + s-ie) 
ln \n + v-s + it)(n + v-s-iey n ' 1 ' 

where a v n) fi v n and 7^ are defined in Eqs.(2.8)-(2.10). By inspection, we see that this 
recurrence relation is deformed to the one in Eq.(2.7) if we systematically redefine the 
coefficients as 

_ n T(v + 1 - s + ie)T{v + 1 - s - ie)T(n + v + 1 + s + ie)T(n + v + 1 + s - ie) 
n ~ 1 Y[v + 1 + s + ie)T(v + 1 + s - ie)T(n + v + 1 - s + ie)T(n + v + 1 - s - ie) ™' 

(3.10) 

where we chose the initial condition 6q = 1. Since the recurrence relation obtained for 
the Coulomb expansion case is identical to the one for the hypergeometric case, the 
renormalized angular momenta v derived from both solutions are the same which allows 
us to relate these two solutions. 

As for the convergence of series in Eq.(3.3), we find 

h u h u eK 

lim = Km n p- = - — , (3.11) 

and from the recurrence relation (3.5) 

lim 1^ {Z ] , = lim = 2 ~, (3.12) 



so that 



lim h nF n+ u{z) = Um b" n F n+v (z) = _e_n 
Thus we find that the series converges for z > ek or | x |> 1. 
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In order to derive the asymptotic behavior of the Coulomb solution R^, it is useful 
to rewrite as 

R C = Rcin + R Cout (3-14) 

where 

R v Cm = e~ lz z u - s (l + ^ ^ 2 u e tw{u+1 ~ s+ie) 

x £ b u n {-2z) n -^\^{n + v+l-s + ie,2n + 2v + 2-2iz). 

n=-oo l(n + z/ + l + s ?ej 

(3.15) 



n Cout — a z 



iz z u-s ^ _|_ 2 ^ T ' ) 2^ e -«7r(^+l+s-ie) 



00 

x £ 6^(-22) n ^(n + z/+ l + s-ie,2n + 2z/ + 2;-2^), (3.16) 

n=— 00 

Another independent solution is obtained by replacing v with —v — \ because —v — \ is 
the renormalized angular momentum if v is the one. Thus, we have another independent 
solution by Rq U ~ 1 ■ The coefficients bz 1 ^ 1 are obtained from b v n by the relation 

bZl- 1 = (-l)X (3-17) 

by choosing 6q = 1 = 1 in conformity with Eq.(3.10). With the use of Eq.(3.14), we 
find by using the identity ^(-L ± s =F ie, -2L, x) = x 2L+1 ^(L + 1 ± s =F ie, 2L + 2, x), 

RcV = -ie- ^ nni r S + i6 \ Rc m , (3-18) 
Cm sunr(v + s - ie) Cm v ' 

Rcout = i^Rcouv (3-19) 
Thus the solution Rq v ~ 1 is expressed by 

r>—v—i — ztti/ gin7r(i/ — s + ie) ->ni/ /o 9n \ 

^ s i n7 r(v + s-te) Rcm + te Rm (3 - 20) 
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4 The relation between two solutions 



First we notice that Rq and R c are solutions of Teukolsky equation. Second we see that 
if we expand these solutions in Laurent series of x — —z/eK, both solutions give the series 
with the same characteristic exponent at x — > oo. Thus, Rq must be proportional to Rq, 

R" = K V R C . (4.1) 

The constant factor K v is determined by comparing like terms of these series. We find 

{eK)- v - r+s 2- v - r (-i) r Y{\ -s-ie-ir) 
T(l + r + v + ir)r(l + r + z/ — s — ie)Y{\ + r + u — s + ie) 



~ r(w + z/ + 1 + ir)r(n + r + 2z/ + 1) „ 
X ^ (n-r)\r(n + is + l-ir) ° n 



n=r 

r ru 



n=— oo 



(r -n)!r(n + r + 2z/ + 2) J ' 

(4.2) 

where r is an arbitrary integer. 

By using these relations, R v in can be written by using the Coulomb expansion solutions 

as 

R in = ( K v R Cin + K-v-lRc'in 1 ) + ( K v R Cout + K -u-lR^out) 

= (K u - ie ^ ^^t±^ K. v . 1 )R^ n + (K u + ij™K_ v _,)R Cout . (4.3) 

The asymptotic behavior at z — > oo is 

iC = AZe-z-^-^ + A^e^z- 1 -, (4.4) 

where A s ^ t and Af£ are amplitudes of the outgoing and incoming waves at infinity of the 
solution which satisfys the incoming boundary condition at the outer horizon. They are 
given by 

oo 

AZt = e-^+^-^-^iKv + ie^K^) £ 6£(-i) n , (4.5) 
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and 



ASU 

A in 



e 2 



7r(— v— ie)Q— ie / ^ 



ze 



-i^ sin7r ( I/ - g + ie ) 

sin7r(z/ + s — ie) 



f ,„. n r(n + i/+l-s + ie) 



(4.6) 



r(n + 1/ + i + s — it) 

One application of these amplitudes is to derive the absorption coefficients. By using 
the method given in Ref.|I4]], the absorption coefficient V can be expressed in terms of 



A™ and A° wt as follows; 



A-sv ASU 
"■out "-out 



A.-SV ASU 

■"■in ■"■in 



(4.7) 



In the end of this section, we show how the upgoing solution which satisfies the 
outgoing boundary condition at infinity is expressed in terms of Rq and Rq v ~ 1 defined 
in Eq.(2.21). From Eq.(2.30), we find 



T>U 
^out 



(A K ir -i^ sin7r ( V ~ S + ie ) A K \r>v 

{A U K V - ie - — ■ —A_ v _ x i^-u-\)ticin 

sin 7T{U + S — 16) 

+ (A U K U + ie™A_ v _ x K_ v _ x )R v Cout . 



By using Eqs.(4.3) and (4.8), we obtain 



T3U 
■^up 



sin ixiv 



ie) 



sin 7r(z/ + s — ie 



-1 D-I/-1 



X 



Jliiru 



+ 



sin ir(v — s + ie) 



Sin 7T v 



ie) 



(4.9) 



5 Low frequency expansions of solutions 

In this section, we discuss how to derive the solution in the expansion of the small 
parameter e = 2Mu. In order to find the solution in Eqs.(2.1) and (4.3) up to some 
power of e, we have to calculate v and a v n to that order by using Eq.(2.12) and (2.13) 
with the condition (2.14) and = a®"' 1 = 1. Other coefficients b v n can be calculated 
from a v n by using the formula (3.10). 
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For a u n with n > 1, the equation for R n (u) is useful. Since a v w 7^ ~ 0(e) and and 
/3£ ~ n(n + 2/ + 1) ~ 0(1), we find R n (u) ~ 0(e) for all positive integer n. As a result 
with dp = 1, we find 

< ~ 0(e") for n > 1. (5.1) 

Before discussing the coeficients for n < 0, we derive the renormalized angular mo- 
mentum v up to 0(e 2 ). For this, it is convenient to use the constraint for n = 1, 
Ri{v)Lq(v) = 1. We notice that Ri(is) ~ 0(e) so that L (u) must behave as 0(l/e). 
which requires that (3q + 7qL_x(z/) ~ 0(e 2 ) because «q ~ 0(e). In order to obtain v 
up to 0(e), we need to know the information of (3q up to 0(e 2 ) where the second order 
term of v involves. Thus, we need the information about R\{y\ L_i(z/), and 7q up 
to 0(e). Here we assume that L_ 2 (^) ~ 0(e) whose validity will be discussed later. In 
this situation, -Ri(^), L_i(u), and 7q can be calculated immediately. By substituting 
these to the constraint eqaution Ri(v)Lq(v) = 1 to find 



v = 1 + 



21 + 1 



s 2 [(Z + l) 2 -s 2 ] 2 (Z 2 -s 2 ) 2 



1(1 + 1) (2/ + l)(2/ + 2)(2/ + 3) (2/ - 1)2/(2/ + 1) 



e 2 + 0(e 3 ). 
(5.2) 



The fact that the correction term of v starts from the second order term of e simplifies 
the calculation of the coeficients up to 0(e 2 ). 

Now we discuss the coeficients for negative integer n for s 7^ which are derived by 
using the equation for L n (v). For large negative value of | n |, L n (v) ~ —itn/2n. Most 
of the negative integer value of n, L n (v) ~ 0(e). There arise some exceptions for certain 
values of n because the denominator of a v n vanishes at n = — I — 1 or — / — | and also (3% 
vanishes at n = —21 — 1 in the zeroth order of e. Because of this, we find for integer Z's, 

L-i-\{v) ~ 0(1) , 

~ 0(l/e) , 

L n {v) ~ 0(e) for all others. (5.3) 

We also find for half-integer n's, 

i- (l+ i H M~o(i/ £ ) , 
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L_ 2{l+h) (v)~0(l/e) 







L n (z/) 


-0(e) 


for all others. 


From the above estimates, 


we find for a integ 


er /, 


< 




O(eH), 


for - 


■ 1 > n > -I, 


a -i-i 


r^j 








< 


rsj 


( e l«l-i) 


, for 


- / - 2 > n > 


V 

a -2l-l 










< 


f-NJ 


0( e l»|-3) 


, for 


- 21 - 2 > n, 


and for a half-integer /, 










< ~ o 


(e |n 


), for 


ra > -(/ 





(5.4) 



(5.5) 



0( e H- 2 ), for - (/ + I) - 2 > „ > -2(1 + ±) - 1 



< ~ 0(e |n| - 4 ), for -2(/ + i)-l>n. (5.6) 

With the above order estimates, we see that how many terms should be needed to cal- 
culate the coeficients with the specified accuracy of e. 

Comming back to u, we assumed that L_ 2 (^) ~ 0(e) whcih is valid if we consider 
/ > §. However, this speciality is due to the fact that we solved the constraint equation 
for n — 1 in Eq.(2.13). Since v is independent of what n we used for solving the constraint 
equation, the result in Eq.(5.2) should be valid for all angular momentum case. In fact 
the result is nonsingular for all integer and half-integer values of /. 

The coefficients a u n and also b v n up to 0(e 2 ) (which are valid for I > |) are obtained 
explicitly by 

.(/ + 1 - s) 2 [(l + l)K + imq}^ 



2(1 + 1) 2 (2/ + 1) 
(I + 1 - s) 2 



2(1 + 1) 2 (2/ + 1) 



(/ + 1)k + imq 2 

1 - t w+-mT2) mq \ 
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e 2 + 0(e 3 ), (5.7) 



„„_ (l + l-s) 2 (l + 2 -s) 2 [(l + 1)k + imq}{(l + 2)K + imq] 2 

4(J + 1) 2 (/ + 2)(2/ + 1)(2/ + 3) 2 6 + ^ j ^ 



. (7 + s) 2 [Ik — imq] 

a - x = 1 2/2(2/ + r 

(Z + S ) 2 



2/ 2 (2/+ 1) 



//t — img 2 

1+ \/-l)/2(/ + l) m9S 



e 2 + 0(e 3 ), (5.9) 



(/-l + s) 2 (/ + s) 2 [(/- l)K-imq][lK-imq] 2 3 
°- 2 = 4(/-l)/ 2 (2/-l) 2 (2/+l) 6 + 0(6 } - (5 - 10) 

The coefficients b v n are given from Eq.(3.10) by 



6 2" (/ + 1 _ s) 2 (/ + 2 _ s) 2 a 2+^(0, (5-12) 

&-i = -^|^^ 1 + 0( e 3 ), (5.13) 

v (i-i-*n-s) 2 a u +0(e3) (514) 

6-2" {l _ l + s) 2 {l + s) 2 a -2 + U ^)- I 5 ' 14 ) 

By using these coefficients, we can evaluate the ingoing and the outgoing amplitudes in 
infinity. From Eq.(4.2), we find by taking r = that K v ~ 0(e~ l+s ). On the other hand, 
the estimate of K_ v _\ needs some care. By taking into account of the singular behaviors 
of gamma functions and the fact that the deviation of v from / starts from the second 
order of e, we find that K_ v _i ~ 0(e l ~ l+s sini^r). Thus we obtain 

~ Oie 21 ' 1 sinmr), (5.15) 

where r = (e — mq)/K. In the approximation up to 0(e 2 ), we can neglect K- v _\ term 
when we restrict / > 3/2. We note that for the Schwarzshild case, r = e so that the ratio 
in Eq.(5.15) is of order e 21 . 
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Thus for Z > 3/2, we get the simple expressions for the outgoing and the incomming 
amplitudes as follows; 

2 



Kit = -<e-w^- fc )2- 1 - +fc ir l/ £ 6£(-i) B , 



(5.16) 



n=-2 



and 



v . n T(n + v + 1 - s + ie) 



(5.17) 



n=-2 T(n + u + l + s-ie)' 

By substituting the coeficients, we can easily calculate the amplitudes up to the order 
e 2 . Since the explicit expressions are complicated, we present the amplitudes up to 0(e) 
explicitly. We find 



and 



^ = _^ e -f^(^-^)2- 1 - s + ie ^ e i[ ^ (1+ ^ )e+02e2l+sh W) e+,/ ' 2e2 



x 



mqs 2 

1 + W(uW e + d2 \ 



(5.18) 



where 



4>2 



A* 



i e ^^(-^+ s -'^)2- 1+s ~ i€ K u e il ~^ {1 ~ T ^ )e+ ' l> ' 2e2] 

mqs 2 , 2 

x 1 + ^ n e + d 2 e 2 



King 



2(2/+ 1) 
1 



(1-8) 



2Z 2 (Z + l) 2 



(Z - 1 - S ) 2 



2(1- l)l 2 (2l- 1) (Z-1)Z 3 (Z + 1) 



+ 



(/_ s )2 ( /+ l + s )2^| f(Z- S )2 (Z+l + s) 1 



^2 = 



4(2/ + 1) [ Z 

+(J->-J-l), 
1 



+ 



Z + l 



2Z + 1 



7 + (mgs)2 I FWTI) + 4P(7Ti? 



i 2 (z + i) 2 

21 + 1 1 



[k(Z- l)Z-(mg) 2 ][(Z-l)Z + s 2 



2(Z-1)Z 2 (2Z-1) 



+ (Z^-Z-1), 



02 = 



mq 
21 + 1 



1 /t5 2 (Z 2 - 5 2 ) K[(Z - l) 2 -5 2 ][Z 2 -5 2 ] 

Z + 2(Z- l)/ 3 (/ + 1) + 4(Z-1)Z 2 (2Z-1) 



1 - 



8Z 2 (Z + 1) 2 [ 1(1 + 1) 



+ (J->-J-l), 



(5.19) 



(5.20) 



(5.21) 



(5.22) 
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(■wigs)' 



Z 2 + s 2 



1 + 



{mqs) 2 
(Z-l)/ 2 (/ + l) 



4/ 2 (Z + l) 2 2/ 2 (2/ + 1) 
[k 2 (1 - 1)1 - (mg) 2 ][((Z - l) 2 + s 2 )(/ 2 + s 2 ) + 4(Z - l)Zs 2 



K 

+ 16 T 



Z(Z + 1) 



4(Z-1)Z 2 (2Z-1) 2 (2/ + 1) 



+ (Z->-Z-l). 



(5.23) 



The above result shows that the absorption coeficient T in Eq.(4.7) is zero up to the order 
of e 2 for Kerr black hole. 



6 Summary and Remarks 

Analytic solutions of Teukolsky equation are obtained in the form of series of hyper- 
geometric functions and Coulomb wave functions. The convergence of these solutions are 
examined. The series solution of hypergeometric type is convergent in the region except 
infinity, while the one of Coulomb type is convergent whem | x |> 1. The renormalized 
angular momentua v turns out to be identical for these two solutions. This fact enabled 
us to relate these two solutions analytically. 

We examined the e dependence of a v n and found that the series corresponds essentially 
to e except for some negative integer n where some anomalous behaviors occured for 
which we need to pay a care to evaluate the coefficients. We explicitly calculated v, the 
coefficients, A s Q u ut and up to the order e 2 . 

The solutions are not only useful to discuss the low fequency behaviors of various 
physical quantities by applying our solutions, but also useful to know the general proper- 
ties of solutions. For example, we consider the e dependence of the renormalized angular 
momentum v in the Schwarzschild metrics, v is determined by solving the transcendental 
equation (2.13) which is composed of f3% and aiklk+i- These quantities are even functions 
of e in the Schwarzschild geometry because r = e. Therefore, we conclude that v is an 
even function of e, i.e., v(—e) = K e )- This property is the special one and not valid for 
the Kerr geometry. The fact that the solutions are given by the e expansion is important 
because the e expansion corresponds to the Post-Minkowskian G expansion and also to 
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the post-Newtonian expansion when they are applied to the gravitational radiation from a 
particle in circular orbit around a black hole. The solutions can be used for the analysis of 
the gravitational radiation from coalescing compact binary systems. Since the analytical 
properties and the convergences are known, the solutions will give a powerful method for 
numerical computation and will contribute to construct the theoretical template towards 
the gravitational observation by LIGO and VIRGO projects. 
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Appendix A: Spheroidal Teukolsky equation 



To expand radial Teukolsky equation, we have to derive the eigenvalue of the spheroidal 
Teukolsky equation. Fortunately, our method is also available in expansion of spheroidal 
Teukolskuy equation. Fackerell expand spheroidal Teukolsky equation in terms of Jacobi 



functions JL2]. In our expansion method, we can derive the eigenvalue of the spheroidal 
field equation which appear in the radial Teukolsky equation. 
The separated spheroidal Teukolsky equation is 

d 2 _ d _ o rn 2 + s 2 — 2msx 



where £ = au, x = cosO. 
We make transformation as 



1 -x 2 



2s£x + E 



S(x) = 0, (A.l) 



S(x) 



1 — x 



1 — X 



U[X) 



(A.2) 



E _ a + (3 ( a + (3 + ^ 



u 



2 ) V 2 

where a = \m + s\, (3 = \m — s\, 
then we recast the equation: 

(1 - x 2 )u" +[(3-a-{2 + a + b)x) u + 

= £ [-2(1 - x 2 )u' + {a + (3 + 2s + 2)xu - (£ + (3 - a)u] . (A.3) 

As a solution of first order equation, we use Jacobi functions which is the solution of the 
equation 

(1 - x 2 )U^ p) " + [(3 -a- (a + (3 + 2)x] U^' + n(n + a + (3 + l)U^ P) = 0. (A.4) 
By using recursion relations of Jacobi functions, we can analytically expand u in terms 

of ui a >®WM, 



u, 



(«./») 

n+j i 



(A.5) 



j=-oo 



where / = n + (a + (3)/ 2. 
We can expand Cj, E in £: 



cj = e c?>e, e = £ E^e, 



(A.6) 



k=0 



k=0 
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here we set c = 1, A = and £ (0) = 1(1 + 1) . 



c , - W 2 -(<*-0) 2 (l c)c (A8) 

C - 1 ~ (2/) 2 (2/+l) (/ (A ' 8) 



and E in Eq.(1.3) which is identical to that of Fackerelljnj who also has shown that the 



convergency of the expansion in terms of Jacobi functions. 

Appendix B: Derivations of equations (2.2) and (3.2) 

The radial Teukolsky equation is written by using the variable y = our with y + = ur + 

and y_ = cur_ as 

d 2 R , 1 1 ,dR 

+ (* + !)(: — — + 



dy 2 y-y+ y-y- dy 

ri l , . t + 2is, l . . e + 2is. 
+ l + (e + is + ) + (e + is ) 

y-y+ « y-y- K 

l (e-is + rf + s 2 l (e - is - r) 2 + s 2 

{y - y+) 2 4 { y - y_y I 

1 e 2 r 2 
+- r(-AH ies - emq)]R = 0. (B.l) 

(y-y + )(y-y-V 22 qn K J 

(a) The expansion in terms of hypergeometric functions 
We define a new variable x by 

V ~ y+ = ck(-x), j/-y_ = e«(l-x). (B.2) 
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We rewrite R as 



R = (-x) a (l - xfR. (B.3) 

in order to eliminate the terms proportional to 1/x 2 and 1/(1 — x) 2 . Then, a and (3 are 
determined to be one of following values, 

a± = ^(-s±i(e-is + T)), (3± = ^(-s ± i(e - is - r)), (B.4) 

respectively. With these choices of ct and /3 and the change of R in the following form 

R = e itKX p, (B.5) 

the equation for p is expressed by 

x(l - x)p" + [(2a + s + 1) - 2(a + (3 + s + l)x]p' 

= —2ieKx(l — x)p + 2ien(a + (3 + 1 + ie)xp + [- 

3 1 

— ien(2a + s + 1) + e«;(e + is) + -e 2 — -r 2 + ies 

where a and 6 are chosen such that the equality 

a + b = 2(a + (3) + 2s + 1 (B.7) 

is satisfied and also they take some simple forms. 

If we take one of choices and the form of R defined in Eq.(B.3) 

becomes a suitable form for the solution which satisfies the incoming boundary condition 
on the outer horizon. On the other hand, the choice of (a + , or (a + , j3 + ) is suitable to 
obtain the solution which satisfies the outgoing boundary condition. In the text, we took 
(ck_,/3 + ) for the solution satisfying the incoming boundary condition in which case the 
above equation (B.6) reduces to the one in Eq.(2.2), by taking a = u+1 — it, b = —u — ir. 
For the solution satisfying the outgoing boundary condition, we took (a + ,f3J) in which 
case the equation becomes a similar form, 
(b) The expansion in terms of Coulomb wave functions 



— abp 

-\ + 2a(3+ (s+ l)(a + (3) 

— emq — ab]p, (B.6) 
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We take the parameterization 

R = [(y- v + )(v - y-)r {1+s)/2 {^ff f ( B - 8 ) 

and determine 7 to eliminate the singularity proportional to l/(y — y-) 2 . Then we find 
7 should take one of the following values 

l± = \[-l±i{e-is-T)\. (B.9) 

If we take one of these values of 7, the equation for / becomes with z = y — y + = —ekx 

z 2 f" + [z 2 + 2(e + i s)z]f = -eKz{f" + /) - 2aenf - ^ - ^ - * S)] f 

z 

+ [A + s(s + 1) - 2e 2 + emq - e/c(e + is)]f. (B.10) 

The equation (3.2) in the text is obtained by taking 7 = 7.. The choice 7 = 7. gives 
a Coulomb type solution which matches with the hypergeometric type solution with 
as we saw in Eq.(4.1). We can also obtain the solution by choosing 7+ which 
matches with the hypergeometric one with 
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